Discontinuities in the level density of small quantum dots under strong magnetic fields 
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Exact diagonalization studies of the level density in a six-electron quantum dot under magnetic 
fields around 7 T ("filling factor" around 1/2) are reported. In any spin-polarization channel, two 
regimes are visible in the dot excitation spectrum: one corresponding to interacting quasiparticles 
(i.e. composite fermions) for excitation energies below 0.4 meV, and a second one for energies 
above 0.4 meV, in which the level density (exponentially) increases at the same rate as in the 
non-interacting composite-fermion model. 
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The lowest energy states of relatively small quan- 
tum dots (number of electrons N e > 3) in strong mag- 
netic fields have been qualitatively described whithin the 
composite-fermion pictureurj. The position of cusps in 
the ground-state energy as a function of the angular mo- 
mentum, and the spin quantum numbers of the low-lying 
excited states are nicely reproduced by this theory. The 
residual interaction between composite fermions (CFs) is 
expected to be a weak, contact interaction. 

The present communication is aimed at giving a quan- 
titative characterization of the density of energy levels 
of small quantum dots in an energy interval below 1 - 
1.5 meV, where dozens or hundreds or states exist. Fully 
converged exact diagonalization results for a 6-electron 
GaAs dot in magnetic fields corresponding to "filling fac- 
tors" v w 1/2 are presented below. The error in comput- 
ing energy eigenvalues is estimated to be lower than 0.02 
meV for levels with excitation energies below 1 meV. The 
studied excitation energy range is small as compared to 
the cyclotronic (Huj c — 12 meV), Coulomb (8.4 meV) or 
confinement {fvjjQ = 3 meV) energies of the model dot. 
Three Landau levels (LLs) are included in the calcula- 
tions, and a 25 meV cutoff in the energy of the non- 
interacting many-electron states used as basis functions 
allows us to deal with hamiltonian matrices of dimen- 
sion less than 850000, which are diagonalized by means 
of a Lanczos algorithm. The main result of the paper 
is the existence of two regimes in the excitation spectra, 
corresponding to low (AE < 0.4 meV) and intermediate 
(AE > 0.4 meV) excitation energies, in which the level 
density increases at different rates. Thus, at an energy 
interval SE around 0.4 meV the level density experiences 
a "discontinuity". 

The model parameters are similar to those used in Ref. 
^. The confinement potential is parabolic. The bare 
Coulomb interaction is weakened by a factor 0.8 to ap- 
proximately account for quasi-bidimensionality (instead 
of exact bidimensionality) . The Zeeman energy is writ- 



ten in the form: 0.01432 B s z meV (B in Teslas and 
s z = ±1/2), corresponding to a 8 nm- width well in mag- 



netic fields around 7 Tea. 
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FIG. 1: The lowest energy levels in each angular momentum 
and spin polarization sectors. Lines are guides to the eyes. 



The lowest energy states in each angular momentum 
and spin polarization tower (the yrast spectrum) for mag- 



2 



netic fields B = 7 and 8 T are shown in Fig. [|. Energy 
jumps between adjacent angular momentum states are 
about 0.6 meV in the spin-polarized case, but roughly 
three times smaller in any other spin polarization sec- 
tor. At this point it is important to stress the role of 
the higher LLs in the energy eigenvalues. The absolute 
contribution of the second and third LLs is around -0.4 
meV, a magnitude much greater than the Zeeman split- 
ting (0.1 meV), and than the characteristic energy spac- 
ing between states near the absolute minimum. Excita- 
tion energies are pushed down 0.1 - 0.3 meV by the higher 
LLs. 
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FIG. 2: The logarithm of the number of levels as a function 
of the excitation energy at B — 7 T: exact results (squares), 
constant temperature fits (dotted lines) , and the shifted NICF 
curves (solid lines). 

We show in Fig. |2| the number of states as a function 
of the excitation energy, AE. In this figure, the reference 
energy is the minimal energy state whithin each polariza- 
tion sector. At low and intermediate excitation energies 
the curves are well fitted by exponential functionsu, 



n(AE) = no exp 



AE 



(i) 



i. e. describedr-by the so called "constant temperature 
approximation"!!! In the polarized case, there are only a 
few states for low excitation energies, thus we do not fit 



them. The fitting parameters no and 6 are given in Table 
|j} Primed quantities correspond to the low-energy sec- 
tors. Notice that the temperature parameter, 0, changes 
very little on going from S z = 2 to S z = for a given 
B. Notice also the apparent jump of the temperature pa- 
rameter at 6E, which means a discontinuity in the level 
density. 
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TABLE I: The level density parameters no and 0. Primed 
quantities correspond to excitation energies below 8E. 

The non-interacting CF (NICF) understanding of the 
excitation spectrum starts from a simplified ILL picture 
in which the excitation energy is written in the form: 



AE = hw r ,„ At 



H u>q/u) c al. 



(2) 



An LL is the variation, with respect to the quasiparticle 
ground state, of the effective LL occupation numbers. 
By definition, AL — \L\ — \ L gs \, where L gs is the ground- 
state angular momentum. fku c is the electron cyclotronic 
energy in GaAs, and the effective cyclotronic energy of 
CFs is extracted from the polarized yrast spectrum, S z = 
3, as: 



E(L 3 + l) + E(L 3 -l)-2E(L 3 ), (3) 



where L3 is the angular momentum of the lowest polar- 



ized state. It gives hui c 



1.15 and 1.19 meV at B = 7 



and 8 T, respectively. The quasiparticles are supossed to 
occupy states in effective Landau levels, which are sepa- 
rated by huj CF . The angular momentum of the quasipar- 
ticle system is L CF which, due to the special form of the 
variational CF wave function, is related to L by 



L= -N e (N e -l) + L c 



(4) 



The NICF model gives qualitatively correct answers 
to questions like the position of cusps in the yrast spec- 
trum, the nature of the first excited states, etc. One may 



3 



ask for other visible manifestations of the NICF in the 
spectrum, and indeed, it gives the correct O parameter 
at intermediate excitation energies, AE > SE. In other 
words, the NICF curve growths exponentially with the 
same 0, although it is shifted with respect to the actual 
spectrum. This fact is illustrated in Fig. || also, where 
the NICF curve at certain level (the sixth in the upper 
figure, for example) is forced to meet the actual curve. 
From this point on the two spectra have the same slope 
in average. 
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the ground and first excited states (low AE), the quasi- 
particles form compact clusterscl and the interaction be- 
tween CFs plays an important role. The low-energy pa- 
rameter O' is a consequence of this interaction. On the 
other hand, at intermediate energies and assuming that 
the interaction is weak, one expects the level density to 
be dominated by the relevant energy scales, i. e. hu; CF 
and h Uq/ll) c . One may visualise this regime in terms of 
Rydberg-like excitations, in which a few non- interacting 
CFs orbit around a core of weakly interacting CFs. En- 
ergy differences will, of course, follow the NICF rules. 

The conclusion to be extracted from this figure is that 
traces of the effective LL structure of CFs may also be 
looked for at intermediate excitation energies, AE > 0.4 
meV. For AE < 0.4 meV the transition to a regime of 
weakly interacting quasiparticles takes place. Tempera- 
ture parameters are different from both sides of SE, thus 
a discontinuity in the level density is expected. 

The level density may be directly measured in Ra- 
man scattering experiments under extreme resonance, 
where the Raman amplitude depends on the density 
of energy levels in final statesB. JVIagnetoconductance 
measurements under equilibrium!!! or non-equilibrium 
conditional, could also give evidence about the level den- 
sity behaviour at intermediate excitation energies. 

Finally, we show in Fig. ^ the magnitudes Tuo CF and 
Os,=3 m a wider magnetic field intensity range, 7 < B < 
12 T. Near B = 12 T, where the "filling factor" is aroud 
1/3, there is a decrease of hw CF which corresponds to 
an abrupt increase of the quasiparticlc massO. A similar 
behaviour is observed in O, showing that the level density 
parameter O depends on a negative power of m CF . 
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FIG. 3: The magnitudes hoJ CF (squares) and Os z =3 (empty 
circles) as a function of magnetic field intensity B. 

There is a natural interpretation of this behaviour. In 
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